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ABSTRACT 



u 

Context. Luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) are key cosmological classes since they account for 
most of the co-moving star formation rate density at z ~ 1 - 2. It is then important to have detailed studies of local samples of their 
j.^ . counterparts for understanding the internal and dynamical processes taking place at high-z. 

Aims. To characterize the two-dimensional morphological, excitation and kinematic properties of LIRGs and ULIRGs we are carrying 
' out an optical integral field spectroscopy (IFS) survey of local (z < 0.26) samples. 

Methods. In this paper we present optical (3800-7200A) IFS with the Potsdam multi-aperture spectrophotometer (PMAS) of the 
northern hemisphere portion of a volume-limited (2750 - 5200 km s~') sample of 11 LIRGs. The PMAS IFS observations typically 
cover the central ~ 5 kpc and are complemented with our own existing HSTfNICMOS images. 
Results. For most LIRGs in our sample, the peaks of the continuum and gas (e.g., Hor, [Nn]/16584) emissions coincide, unlike what 
is observed in local, strongly interacting ULIRGs. The only exceptions are galaxies with circumnuclear rings of star formation where 
the most luminous Ha emitting regions are found in the rings rather than in the nuclei of the galaxies, and the displacements are 
well understood in terms of differences in the stellar populations. A large fraction of the nuclei of these LIRGs are classified as 
LINER and intermediate LINER/HII, or composite objects, which is a combination of starformation and AGN activity. The excitation 
^ — , conditions of the integrated emission depend on the relative contributions of Hn regions and the diffuse emission to the line emission 

over the PMAS FoV. Galaxies dominated by high surface-brightness H n regions show integrated H n-like excitation. A few galaxies 
show slightly larger integrated [Nii]/16584/Ha' and [S n]/16717,6731/Ha line ratios than the nuclear ones, probably because of more 
contribution from the diffuse emission. The Ha velocity fields over the central few kpc are generally consistent, at least to first order, 
with rotational motions. The velocity fields of most LIRGs are similar to those of disk galaxies, in contrast to the highly perturbed 
fields of most local, strongly interacting ULIRGs. The peak of the Ho- velocity dispersion coincides with the position of the nucleus 
and is likely to be tracing mass. All these results are similar to the properties of z ~ 1 LIRGs, and they highlight the importance of 
detailed studies of flux-limited samples of local LIRGs. 
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1. Introduction 

In recent years deep cosmological surveys have been extremely 
successful in identifying large samples of high-z galaxies using 
specific wavelengths (e.g., UV, optical, infrared, submillimeter) 
or combinations of them. Luminous infrared galaxies (LIRGs, 
with infrared 8-1000 yt/m luminosities Lm — 10^ — 10'^L©, 
see Sanders & Mirabel 1996) and ultraluminous infrared galax- 
ies (ULIRGs, with IR luminosities L IR = 10 12 - 10 13 L©, see 
Lonsdale, Farrah, & Smith 2006) are significant cosmological 



* Based on observations collected at the German-Spanish 
Astronomical Center, Calar Alto, jointly operated by the Max- 
Planck-Institut fur Astronomie Heidelberg and the Instituto de 
Astrofisica de Andalucfa (CSIC). 



classes. These IR-selected galaxies are the main contributors to 
the co-moving star formation rate density of the universe at z > 1 
(Elbaz et al. 2002; Le Floc'h et al. 2005; Perez-Gonzalez et al. 
2005; Caputi et al. 2007). 

Because these high-z samples are so numerous, most efforts 
concentrate on characterizing their integrated properties, such 
as stellar masses, star formation rates, average ages, and metal- 
licities. However, to understand fully how galaxies formed and 
evolved, one needs spatially resolved information to study their 
kinematics, stellar populations, the star, gas, and dust distribu- 
tions, as well as the gas excitation conditions. At high-z this has 
been done for optically /UV selected galaxies using integral field 
spectroscopy (IFS). These works find clumpy Ho- morpholo- 
gies with relatively well-ordered velocity fields that are consis- 
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tent with the presence of large disks at z ~ 2, while other sys- 
tems are better explained as merger candidates (see e.g., Forster- 
Schreiber et al. 2006; Genzel et al. 2006, 2008; Wright et al. 
2009). A similar result is found at intermediate redshifts (Puech 
et al. 2008; Yan et al. 2008). 

In the local universe LIRGs and ULIRGs are much less nu- 
merous than at high-z, and a large amount of work has already 
been done to characterize their properties using optical long-slit 
spectroscopy (e.g., Heckman, Armus, & Miley 1987; Armus, 
Heckman, & Miley 1989; Kim et al. 1995, 1999; Veilleux et 
al. 1995, 1999; Wu et al. 1998; Heckman et al. 2000; Rupke, 
Veilleux & Sanders 2005; Chen et al. 2009). The majority of 
optical and near-IR IFS works have so far focused on small sam- 
ples of (U)LIRGs or individual galaxies (e.g., Arribas, Colina, 
& Clements 2001; Murphy et al. 2001; Lfpari et al. 2004a,b; 
Colina, Arribas, & Monreal-Ibero 2005; Monreal-Ibero, Arribas, 
& Colina 2006; Garcfa-Marin et al. 2006; Reunanen, Tacconi- 
Garman, & Ivanov 2007). The work of Shapiro et al. (2008) has 
recently highlighted the importance of having local templates 
for understanding the internal and dynamical processes taking 
place at high-z. We note, however, that care is needed when 
comparing local IR luminosity matched galaxies with high-z sys- 
tems. For instance, at z ~ 2 the mid-IR spectra of star-forming 
ULIRGs are more similar to those of local starbursts and LIRGs 
than to those of local ULIRGs (Farrah et al. 2008; Rigby et al. 
2008; Alonso-Herrero et al. 2009). One possible explanation is 
that in high-z ULIRGs star formation is taking place over larger 
scales, a few kpc, than in local ULIRGs where most of the in- 
frared emission arises from sub-kpc scale regions (e.g., Soifer 
et al. 2000). In contrast intermediate redshift (z ~ 0.7) LIRGs 
appear to be experiencing similar starburst phases to those of lo- 
cal LIRGs (Marcillac et al. 2006). We thus need to understand 
the spatially-resolved physical processes and properties of local 
(U)LIRGs and compare them with those of distant IR-selected 
galaxies. 

We recently started an optical IFS survey of a representa- 
tive sample of approximately 70 nearby (z < 0.26) LIRGs and 
ULIRGs (see Arribas et al. 2008). The general goal of this sur- 
vey is to provide a two-dimensional, as opposed to integrated 
or nuclear, characterization of the physical and dynamical pro- 
cesses taking place in local LIRGs and ULIRGs. As discussed 
in detail by Arribas et al. (2008), we are conducting this sur- 
vey using three different optical IFS instruments in both the 
northern and the southern hemispheres. These are: VIMOS (Le 
Fevre et al. 2003) on the VLT, the INTEGRAL+WYFFOS sys- 
tem (Arribas et al. 1998; Bingham et al. 1994) on the William 
Herschel Telescope, and the Potsdam multi-aperture spectropho- 
tometer (PMAS, Roth et al. 2005) instrument on the 3.5 m tele- 
scope at Calar Alto (Spain). In this paper we present an atlas 
of the IFS observations obtained with PMAS for the northern 
hemisphere portion of the volume-limited sample of LIRGs of 
Alonso-Herrero et al. (2006), which is part of the larger survey 
of Arribas et al. (2008). Additionally, we are observing a rep- 
resentative sample with near-infrared IFS (Bedregal et al. 2009) 
using SINFONI (Eisenhauer et al. 2003) on the VLT. 

The paper is organized as follows. Sect. 2 gives details on the 
sample, the observations and the data reduction. Sect. 3 presents 
the analysis of the PMAS IFS data. Sects. 4, 5, and 6, discuss 
the general results on the morphologies of the emission lines 
and continuum, the nuclear and integrated ID spectra, and the 
kinematics of the ionized gas, respectively. The discussion and 
summary of this work are given in Sect. 7. Throughout this paper 
we used Hq = 75kms~' Mpc~'. 



2. Observations 

2.1. The Sample 

We observed the majority of the northern hemisphere galaxies 
(see Table 1) from the volume-limited representative sample of 
local LIRGs of Alonso-Herrero et al. (2006). This sample was 
drawn from the IRAS Revised Bright Galaxy Sample (RBGS, 
Sanders et al. 2003). The Alonso-Herrero et al. (2006) sample 
is limited in distance (velocities of 2750 - 5200 km s _1 or dis- 
tances of d ~ 35 - 75 Mpc for the assumed cosmology) so 
that the Paa emission line at rest-frame wavelength 1.875/im 
could be observed with the NICMOS F190N narrow band fil- 
ter (see Sect. 2.3 and Alonso-Herrero et al. 2000a, 2006 for full 
details). Arp 299, the most luminous system in this LIRG sam- 
ple, was observed with the IFS INTEGRAL instrument and an- 
alyzed in detail by Garcia-Marm et al. (2006). For the sake of 
completeness we will include Arp 299 (log(LiR/L ) = 11 .88) in 
the discussions presented in Sects. 4, 5, and 6. Our sample of 
LIRGs covers a range in infrared luminosities of log(Li R /L ) = 
1 1 .05 - 1 1 .88. The average infrared luminosity of the full sample 
(northern and southern hemispheres) is Lir = 2.1 x 1O U L or 
log(L IR /L ) = 11.32. Except for Arp 299, which is a strongly 
interacting system, the rest of the northern hemisphere LIRGs 
in this sample are apparently isolated galaxies, weakly inter- 
acting systems (e.g., NGC 7469, NGC 7771) and galaxies with 
small companions but no clear morphological signs of interac- 
tion (e.g., NGC 6701, see Marquez et al. 1996). 

2.2. PMAS Optical Integral Field Spectroscopy Observations 

We obtained optical IFS of 1 1 LIRGs using PMAS on the 3.5 m 
telescope at the German-Spanish Observatory of Calar Alto 
(Spain) during three observing runs: November 2005, May 2006, 
and December 2006. The PMAS observations were taken with 
the Lens Array Mode configuration which is made of a 16 x 16 
array of microlenses coupled with fibers called hereafter spax- 
els. We used the 1" magnification which provides a field of view 
(FoV) of 16" x 16". We used the V300 grating with a dispersion 
of 1.67 A pixeL 1 and an approximate spectral range of 3400 A. 
The wavelength range covered by the observations was approx- 
imately 3800 - 7200 A. The approximate spectral resolution of 
the spectra in the 2x2 binned mode is 6.8 A full width half 
maximum (FWHM, see also Sect. 3.1). 

2.2.1. Observing procedure 

The total integration time for each galaxy was split into three or 
four individual galaxy exposures (see Table 1), each of which 
having between 400 and 1000 s. Given the relatively small FoV 
of PMAS in the Lens Array configuration and the large extent 
of the galaxies, we obtained a separate sky integration for each 
galaxy, interleaved with the galaxy observations and with a com- 
parable single exposure time of 900 or 1000 s. For all the galax- 
ies we obtained one pointing, except for NGC 777 1 for which we 
took two pointings to cover the approximate central 28" x 16" 
region. 

Table 1 gives for each target the observing campaign, inte- 
gration times for the individual galaxy exposures and number of 
exposures, integration time, seeing conditions for each data set, 
as well as the airmasses of the observations. The typical seeing 
conditions of the observations varied between 0.8" and 1.7", 
depending on the observing campaign. 



Alonso-Herrero et al.: PMAS integral field spectroscopy of local LIRGs 
Table 1. Log of the PMAS Observations. 



3 



Galaxy 


Dist 


log Lir 


Date 


Galaxy t m 


Sky t mt 


Airmass 


Seeing 


Conditions 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 23 


59.6 


11.05 


7 Nov 2005 


3 x 1000 


1000 


1.06 


1.0- 1.2 


P 


MCG +12-02-001 


64.3 


11.44 


18 Dec 2006 


4x800 


1000 


1.47 


1.7 


NP 


UGC 1845 


62.0 


11.07 


8 Nov 2005 


3x900 


900 


1.05 


1.4-2.0 


NP 


NGC 2388 


57.8 


11.23: 


7 Nov 2005 


3 x 1000 


1000 


1.04 


1.3 - 1.1 


P 


MCG +02-20-003 


67.6 


11.08: 


17 Dec 2006 


4x 800 


1000 


1.60 


1.2 


NP 


IC 860 


59.1 


11.17: 


29 May 2006 


4x 800 


1000 


1.03 


1.7 


NP 


NGC 5936 


60.8 


11.07 


30 May 2006 


4x800 


1000 


1.17 


1.5 


NP 


NGC 6701 


56.6 


11.05 


30 May 2006 


1 x 400 + 4 x 600 


800 


1.09 


1.3 


NP 


NGC 7469 


65.2 


11.59 


8 Nov 2005 


3x900 


900 


1.15 


1.2 


NP 


NGC 7591 


65.5 


11.05 


8 Nov 2005 


3 x 1000 


1000 


1.32 


1.1 - 1.2 


NP 


NGC 777 1-E 


57.1 


11.34 


7 Nov 2005 


2 x 500 + 2 x 1000 


1000 


1.15 


0.8 


P 


NGC 777 1-W 






7 Nov 2005 


3 x 1000 


1000 


1.06 


0.8 


P 



Notes. — Column (1): Galaxy. For NGC 7771 we took two different pointings to cover the approximate central 28" x 16" region of the galaxy, 
as explained in the text. Columns (2) and (3): Distance (in Mpc) and IR luminosity (in L Q ) taken from Sanders et al. (2003) for the IRAS 
RBGS. Column (4): Date of the observations. Column (5): Integration time (in seconds) of the galaxy observations after discarding bad data sets. 
Column (6): Integration time (in seconds) of the sky observations. Column (7): Airmass at the beginning of the observation. Column (8): Typical 
seeing conditions (in arcseconds) of the observations. If two values of the seeing are given they correspond to the beginning and the end of the 
observations. Column (9): Observing conditions. P: Photometric. NP: Non-photometric. 
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Fig. 2. PMAS maps of the Hj3 and Ha emission lines of 
NGC 7469 fitted using two components, broad (left panels) and 
narrow (right panels). Orientation is north up, east to the left. 
The images are shown on a square root scale. 



During each night we obtained calibration arcs and lamps, 
and fiat-fields. Since PMAS is located at the Cassegrain focus, it 
is affected by the changing flexures of the telescope when track- 
ing the targets. Therefore, we obtained individual arc and inter- 
nal lamps exposures for each different pointing, when applicable 
before and after culmination of the target. The lamps used were 



a ThAr lamp for the Nov 2005 and May 2006 campaigns, and a 
HgNe lamp for the Dec 2006 campaign. Additionally through- 
out the nights we observed spectrophotometric standard stars to 
correct for the instrument response, and to flux calibrate the data. 
We note, however that conditions were non-photometric for sev- 
eral nights in all the campaigns (see Table 1). 

2.2.2. Data reduction 

The PMAS data were reduced using a set of customized scripts 
under the IRAF0 environment. The first step of the data reduc- 
tion was to determine and subtract the bias level. Next, using 
reference internal continuum lamp exposures, we identified and 
traced the location of each of the 256 spectra along the disper- 
sion direction of the CCD. Once we extracted the individual 
spectra, the third reduction step was the wavelength calibration, 
which we carried out using a model obtained from arc lamp ex- 
posures. We checked the wavelength calibration against known 
sky emission lines and measured a median standard deviation of 
1.8, 2.0, and 1.8A for the Nov 2005, May 2006, and Dec 2006 
observing runs, respectively. The fourth step was to correct for 
the sensitivity variations by creating a response image using in- 
ternal calibration lamp images and a sky flat exposure. After that, 
a relative flux calibration was performed using standard star ob- 
servations. The next step was to combine the different galaxy 
exposures (a minimum of three) for each individual pointing to 
improve the signal-to-noise ratio and to remove cosmic rays. The 
sky subtraction was done on a spaxel-by-spaxel basis using for 
each galaxy its own sky observation. The final step was to build 
the data cubes and rotate them to the north up, east to the left 
orientation. 

2.3. HST/NICMOS Observations 

The ffiiT/NICMOS observations were obtained with the 
NICMOS NIC2 camera, which has a pixel size of 0.076" and a 



1 IRAF software is distributed by the National Optical Astronomy 
Observatory (NOAO), which is operated by the Association of 
Universities for Research in Astronomy (AURA), Inc., in cooperation 
with the National Science Foundation. 
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Table 2. Nuclear and Integrated observed (not corrected for extinction) line ratios. 



Galaxy 


Type 


Size 


[oni]/H/; 


[OI]/Ha 


[NIIl/Ho- 


[SII]/Ha 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC 23 




98 


1 59 






n 48 




Integrated 


4.5 


0.58 


0.044 


0.57 


0.33 




M 1 1 1 & q r 
1>I L1L 1 L il 1 


U.J 1 


n (\i 

u.u / 


f) f)99 


n aa 


U.JU 




Integrated 


5.8 


0.61 


0.030 


0.42 


0.35 


unr 1 845 


Ml If If 1 Ql" 




2.07 




1 09 

1 .U7 






Integrated 


3.2 


1.33 


0.065 


0.72 


0.28 


NGC 2388 


Nuclear 


0.28 


0.32 


0.029 


0.63 


0.22 




Integrated 


5.1 


0.60 


0.066 


0.56 


0.30 


MCG +02-20-003 


Nuclear 


0.33 


0.67 


0.054 


0.63 


0.37 




Integrated 


6.5 


0.65 


0.046 


0.45 


0.43 


IC 860 


Nuclear 


0.29 






7.81 


3.85 


NGC 5936 


1 > LK Ik. tl 1 


0.30 


0.32 


0.031 


0.63 


0.22 




Integrated 


4.7 


0.24 


0.034 


0.48 


0.25 


Mf T r fS701 

1 >l Vlv . U/Ul 


Mlir , 1f 1 Q1" 


0.27 


u.uo 


n nfi5 


79 


35 




Integrated 


5.2 


0.62 


0.054 


0.67 


0.35 


NGC 7469* 


Nuclear 


0.32 


9.96 


0.13 


4.10 


0.39 




Integrated 


6.1 


3.53 


0.048 


0.55 


0.30 


NGC 7591 


Nuclear 


0.32 


0.98 


0.11 


0.96 


0.41 




Integrated 


5.1 


0.85 




0.85 


0.59 


NGC 7771 


Nuclear 


0.28 


0.22 


0.032 


0.40 


0.21 




Integrated 


8.0 


0.42 




0.55 


0.43 



Notes. — Column (2): PMAS type of ID spectra. Column (3): Linear physical size in kpc covered by the PMAS extraction apertures for the 
nuclear and integrated spectra using the distances given by Sanders et al. (2003) for the IRAS RBGS. Columns (4), (5), (6), and (7): Observed (not 
corrected for extinction and/or stellar absorption) line ratios. 'For NGC 7469 the line ratios correspond to the narrow components of the hydrogen 
recombination lines. 



FoV of ~ 19.5" x 19.5". Details on the observations and data re- 
duction procedures can be found in Alonso-Herrero et al. (2006). 
For this paper we make use of the 1 .6 fim continuum observa- 
tions and the continuum-subtracted Paa images for comparison 
with the PMAS data. The only additional step needed for the 
NICMOS images was to rotate and trim them to match the ori- 
entation and FoV, respectively, of the PMAS images. The angu- 
lar resolution of the NICMOS data is approximately 0.15" and 
0.18" (FWHM) for the 1.6/im continuum and Paa images, re- 
spectively. 

3. IFS Data Analysis 

3. 1 . Spectral maps 

We constructed spectral maps of the brightest emission lines 
by fitting the lines to Gaussian functions and the adjacent con- 
tinuum to straight lines, on a spaxel-by-spaxel basis. To do so 
we developed our own routines which make use of the IDL- 
based MPFITEXPR algorithm^ developed by Markwardt (2008) 
to measure in an automated fashion the central wavelength, the 
width of the Gaussian (cr = FWHM/2.35), and integrated flux 
for each emission line. Using this algorithm we were also able to 
fix the relative wavelengths and line ratios according to atomic 

2 http://www.purl.com/net/mpfit 



parameters when fitting multiple emission lines ([Nn]/16548, 
Ha, [Nn]^6584, and [S n]/167 17,6731). Additionally for each 
of the two sets of lines we imposed that the lines had the same 
width. We did not attempt to correct the maps of the H/3 emission 
line for the presence of H/i in absorption, which is observed in 
most of the spectra of the sample of LIRGs (see also Sect. 3.2). 

The spectral maps of the emission line fluxes, as well as the 
maps of the central wavelength (or velocity field), and cr (or ve- 
locity dispersion) were only constructed for those spaxels whose 
integrated line flux was 3 6 above the local continuum, where 
5 was the standard deviation of the fitted local continuum. We 
found that the emission lines of all galaxies were adequately fit- 
ted with one component (Fig. 1), except for NGC 7469, a well- 
known Seyfert 1 galaxy, for which we fitted the hydrogen re- 
combination lines using a broad and a narrow component. The 
maps of both components are shown in Fig. 2. 

The uncertainties of the measured velocities and velocity 
dispersions depend on the errors in measuring the centroid and 
width of the emission line, which in turn depend on the signal- 
to-noise of the spaxel, and the systematic errors (e.g., the wave- 
length calibration, see Sect. 2.2). The typical errors in measur- 
ing the widths of the Gaussians when fitting the Ha+[Nn] lines 
were 10%. The maps of the observed Ha velocity dispersion 
(o" bs) were corrected for the instrumental resolution (cr; ns ) us- 
ing cr 2 = cr 2 bs - cr 2 ns . For each observing campaign we measured 
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Fig. 3. (a) Maps of the observed Ha velocity field (left panels) and the velocity dispersion (right panels). The latter has been corrected 
for the instrumental resolution. The zero points for the Ha velocity fields are set at the location of the peak of the 6200 A continuum 
emission as marked with the crosses for each galaxy. Orientation is north up, east to the left. 



the instrumental resolution on a spaxel-by-spaxel basis using the 
observations of the arc emission lines at a wavelength close to 
that of Ha. The median instrumental resolution across the detec- 
tor near Ha is cr ins = 2.9+0.2 A for all three observing runs. The 
Ha velocity fields and the maps of the Ha velocity dispersion are 
shown in the left and right panels, respectively, of Fig. 3. 

Since the emission lines observed in IC 860 are faint, we 
decided to fit the Ha+[Nn] and the [S n] lines (H/3, [Om]^5007 
and [Oi]/I6300 were not detected) manually using splot within 
iraf without any constraints. The spectral maps for this galaxy 
are shown in Fig. If. As can be seen from this figure, the Ha and 
[S n] lines are only detected in the nuclear regions, whereas the 
[Nn]/16584 emission line is more extended. 



In addition to fitting the brightest emission lines for each 
galaxy, we constructed a continuum image centered at A st 
6200 A by summing up the continuum spectra over a rectangular 
band width of approximately 17 A. Fig. 1 shows the continuum 
maps for our sample of LIRGs. The peak of the optical con- 
tinuum emission is marked with a cross for each galaxy on the 
PMAS maps in this figure as well as in Figs. 2 and 3. 

As explained in Sect. 2.1.1, for NGC 7771 we took two 
separate pointings to cover the approximate central 28" x 16" 
region of the galaxy. The systematic error of telescope offsets 
with PMAS is given by the astrometric accuracy of the PMAS 
Acquisition and Guiding (A&G) system, which has been deter- 
mined to be 0.1966 ± 0.0004"/pixel (Roth et al. 2005). As the 
A&G system is an integral part of the instrument, systematic 
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offset errors are negligible. The statistical accuracy of the auto- 
guider is a fraction of a pixel, and the canonical value is 0.1" 
r.m.s. (Roth, 2009, private communication). Thus, we used the 
offsets commanded using the PMAS acquisition images to con- 
struct mosaics of the emission line and continuum maps, as well 
as the Ha velocity fields and maps of the velocity dispersion. 

Fig. 1 shows for the 1 1 LIRGs in our sample the PMAS maps 
of the brightest optical emission lines, the HS77NICMOS con- 
tinuum subtracted Paa images, together with images of the stel- 
lar emission at ~ 6200 A (PMAS) and at 1.6 //m (NICMOS). 
Since the PMAS images do not have astrometry, the PMAS and 
NICMOS images were registered by using the peaks and shapes 
of the continua for reference. For the typical distances of our 
galaxies the FoV of the PMAS observations cover the central 
4.3 -5.3 kpc, except for NGC 7771 for which the PMAS mosaics 
cover approximately the central 7.8 kpc x 4.5 kpc. The maps are 



shown on a square root scale to maximize the contrast between 
diffuse and bright regions. 



The observed Ha velocity fields for all the LIRGs in our 
sample except for IC 860 are shown in Fig. 3. For NGC 7469 
we show the [Nn]/16584 velocity field instead, as that of Ha 
is affected by the uncertainties associated with fitting the broad 
and narrow components. The zero points of the velocity fields 
are set at the peak of the 6200 A continuum emission. The only 
exception is NGC 7771 where the zero point is placed at the po- 
sition that makes the velocity field gradient symmetric, and it is 
approximately coincident with the peak of the near-infrared con- 
tinuum. Table 3 gives the measured Ha cz for the nuclei of our 
sample of LIRGs. 
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Fig. 3. (b) As Fig. 3a. The velocity field and the map of the velocity dispersion for NGC 7469 are for the [Nn]/16584 emission line 
instead of Ha. 



3.2. Extraction of 1 D spectra 

For each galaxy we extracted two ID spectra: the nuclear spec- 
trum and the integrated spectrum. The nuclear spectra corre- 
spond to the spaxel at the peak of the 6200 A continuum emis- 
sion. The physical size covered by the nuclear spectrum is given 
for each galaxy in Table 2, and it is typically the approximate 
central 300 pc. For reference the nuclear spectra of Kim et al. 
(1995) and Veilleux et al. (1995) were extracted with a linear 
physical size of 2 kpc (see discussion in Sect. 5). 

The integrated spectrum of each galaxy was extracted by 
defining ~ 6200 A continuum isophotes and then summing 
up all the spaxels contained within the chosen external contin- 
uum isophote. The external isophotes (plotted for all galaxies in 
Fig. 1) were selected to cover the PMAS FoV as much as pos- 
sible, without compromising the quality of the extracted spec- 



tra. In Table 2 we give for each galaxy the approximate phys- 
ical size along the major axis of the galaxy of the outer con- 
tinuum isophote used for the integrated spectra. Note that the 
term integrated is used in the sense of integrated spectra over 
the PMAS FoV, and the integrated spectra do not encompass the 
whole galaxy (see e.g. figure. set. 8 of Moustakas & Kennicutt 
2006 for a comparison with one of the galaxies in our sample, 
NGC 23). Typically our integrated spectra cover the central 3 
to 8 kpc along the major axis of the galaxy, depending on the 
galaxy, but for most galaxies they include the central ~ 5 kpc 
(see Table 2). The integrated spectra of those galaxies observed 
under non-photometric conditions (Table 1 ) may be affected by 
imperfect sky subtraction, but this does not affect the measure- 
ments of the brightest emission lines. The full nuclear spectra are 
shown in Fig. 4 for each LIRG in the sample. In the same figure 
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we present in the insets the blue part of the integrated spectra 
to emphasize the differences between the underlying absorption 
features and the H/3 and [Oiii]/15007 emission lines in the nu- 
clear and integrated spectra. 

The fluxes of the brightest emission lines of the nuclear and 
integrated spectra were measured manually using splot within 
iraf and were compared with the automated flux measurements 
used for constructing the spectral maps (see Sect. 3.1), except 
for NGC 7469. In the case of the manual measurements we 
did not impose any constraints on the line widths and ratios 
when measuring the [Nn]/16548, Ha and [Nn]/16584, and the 
[S n]/Ll6717,6731 lines. For the nuclear and integrated spectra 
of NGC 7469 we used the method and restrictions described in 
Sect. 3.1. We find a good agreement between the manual and 
the automated measurements of the nuclear values. The largest 
differences (up to ~ 30%) are for the [Oi]/I6300/Ha line ratio, 
whereas for the other line ratios the differences are always of less 
than 15%. As done for the spectral maps, we did not attempt to 
correct for the presence of H/3 in absorption. The observed (not 
corrected for extinction) line ratios for the nuclear and integrated 
spectra are given in Table 2. 

4. Morphology of the stellar and gas emissions 

The optical and near-infrared continuum images (Fig. 1) reveal 
the presence of bright nuclei, and a large number of star clus- 
ters in the nuclear regions as well as along the large scale spiral 
arms. The star clusters are only unveiled by the higher angular 
resolution of the NICMOS images, which is typically a few tens 
of parsecs for our sample of LIRGs. There is a good overall cor- 
respondence, on scales of a few hundred parsecs (the PMAS spa- 
tial resolution), between the optical and the near-infrared stellar 
continua as mapped out by the PMAS 6200 A and the NICMOS 
1 .6 /J.m emissions, respectively. This indicates that for the ma- 
jority of these LIRGs the effects of extinction on the contin- 
uum morphologies are not overly severe, except for the nuclei 
of the galaxies (Alonso-Herrero et al. 2006). This is clearly seen 
in NGC 7771 where the peaks of the optical and near-infrared 
continua appear displaced by a few arcseconds (see Figs. Ik 
and 3). This is probably due to the highly inclined nature of this 
galaxy, as well as the diversity of stellar populations and patchy 
extinction present in the ring of star formation (Davies, Alonso- 
Herrero & Ward 1997; Smith et al. 1999; Reunanen et al. 2000). 
Another example with large differences between the optical and 
near-infrared continuum emission is Arp 299, and we refer the 
reader to Alonso-Herrero et al. (2000a) and Garcia-Marm et al. 
(2006) for a full discussion. 

The PMAS Ha and the NICMOS Paa emissions are well 
correlated, and both trace the nuclear emission as well as 
the emission from bright, high surface-brightness H n regions 
(Fig. 1). The high angular resolution of the NICMOS Paa im- 
ages resolves with exquisite detail the sites of the youngest star 
forming regions in the central regions of LIRGs. These high 
surface-brightness H n regions can be either located in the cen- 
tral 1-2 kpc or spread out throughout the disk of the galaxies 
(see Alonso-Herrero et al. 2006 for more details). 

Since the NICMOS Paa and nearby continuum images were 
taken with narrow-band filters, and a relatively small pixel size, 
they are not very sensitive to the diffuse low surface-brightness 
emission. This may also be due to the fact that in LIRGs the dif- 
fuse emission suffers much less extinction than the bright H n re- 
gions (see Rieke et al. 2009 and references therein). The PMAS 
Ha images show emission from the H n regions, at lower angu- 
lar resolution than the NICMOS images. Additionally, for most 



LIRGs the PMAS Ha images are also sensitive to more dif- 
fuse low surface-brightness emission. This extended emission 
can be seen almost over the entire PMAS FoV (e.g., NGC 5936 
Fig. lg and NGC 6701 Fig. lh), and beyond, as shown by other 
works for a few galaxies in common with our sample (Marquez, 
Masegosa, & Moles 1999; Dopita et al. 2002; Hattori et al. 
2004). It is only the two galaxies with the most compact nuclear 
Paa emission, IC 860 (Fig. If) and to a lesser degree UGC 1845 
(Fig. lc), that also show relatively compact Ha emission. 

The fact that Ha and Paa morphologies are in general similar 
suggests that over the PMAS FoV and with the PMAS angular 
resolution the extinction effects on Ha are not severe, except in 
the very nuclear regions. Indeed, Alonso-Herrero et al. (2006) 
for the galaxies in common with this work measured average 
extinctions to the gas of Ay ~ 2 - 4 mag over the Paa emit- 
ting regions (a few kpc). Similar values of the extinction were 
found by Veilleux et al. (1995) from the Balmer decrement. The 
PMAS Hy6 maps are similar to those of Ha although they are 
more affected by extinction and/or the presence of H/3 in absorp- 
tion, especially in the nuclear regions (e.g., UGC 1845 Fig. lc, 
MCG+02-20-003 Fig. le, NGC 7591 Fig. lj). 

In general the overall morphology of the brightest forbidden 
lines ([Nn]/16584, [S ii]/1/16717,6731) shows a reasonable corre- 
lation with that of Ha. On smaller scales, however, some differ- 
ences are already apparent in Fig. 1. For instance, in NGC 23 
(Fig. la) all the optical emission lines except for Ha and H/3 
peak in the nucleus, whereas the brightest regions of hydrogen 
recombination line emission (Ha, Paa), which trace the sites of 
on-going star formation, are in the circumnuclear ring of star for- 
mation. There are also small scale differences between the Ha 
and [Nn]/16584 emissions of the nuclei and the Hn regions For 
example, some of the extra-nuclear Hn regions of NGC 2388 
(Fig. Id) and MCG +02-20-003 (Fig. le) appear to have lower 
[N n]/l6584/Ha ratios than their nuclei. This is a well known 
behavior of galaxies (Veilleux & Osterbrock 1987; Kennicutt, 
Keel, & Blaha 1989; Sarzi et al. 2007). The small scale differ- 
ences of the different emission lines in LIRGs are more apparent 
from the spatially resolved properties of the optical line ratios of 
this sample of LIRGs and will be discussed in detail in a forth- 
coming paper (Alonso-Herrero et al. 2009, in preparation). 

Although the limited angular resolution of the PMAS contin- 
uum and gas maps does not allow us to resolve all the morpho- 
logical details seen in the NICMOS images, the general behav- 
ior is the same. For the majority of the LIRGs in our sample the 
peaks of the stellar and the Ha gas emissions are coincident and 
located in the nuclear region. The few cases of displacements be- 
tween the stellar and the Ha peaks are found in those LIRGs with 
circumnuclear rings of star formation without Seyfert activity. In 
NGC 23 (Fig. la) and NGC 7771 (Fig. Ik) the stellar emission 
peaks in the nuclei, whereas the brightest Ha and Paa emission 
are in luminous Hn regions in the rings. These displacements 
are likely to be due to differences in the stellar populations. That 
is, the youngest regions and thus brightest Ha emitting regions 
are in the rings, whereas the nuclei contain older stellar pop- 
ulations. The latter is clearly demonstrated by the presence of 
strong absorption features (Balmer line series, H+K Call lines) 
in the nuclear spectra of these two galaxies (Fig. 5). In the case 
of NGC 7469, which also shows a bright circumnuclear ring of 
star formation (see Genzel et al. 1995; Diaz-Santos et al. 2007, 
and references therein), the peaks of the stellar and gas emission 
are coincident with the nuclear AGN. 

In contrast to the majority of the LIRGs in this volume- 
limited sample, in the interacting LIRG Arp 299 the peaks of the 
observed warm ionized gas are displaced from the stellar peaks 
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Fig. 3. (c) As Fig. 3a. In the case of NGC 7771 the asterisk marks the position that makes the velocity field gradient symmetric, and 
it is approximately coincident with the peak of the near-infrared continuum and the maximum of the velocity dispersion. 



typically by ~ 1 .4 kpc (see Garcfa-Marfn et al. 2006). In local in- 
teracting ULIRGs the displacements between the continuum and 
gas emissions are common and even larger, typically 2-4 kpc 
and in some exceptional cases as large as ~ 8 kpc (Colina et al. 
2005; Garcfa-Marfn et al. 2009). The differences in the stellar 
and gas distributions in ULIRGs are understood in terms of the 
more extreme effects produced by the interaction processes on 
the spatial distribution of the ionizing sources and the presence 
of large amounts of dust in the nuclear regions. 



5. Nuclear versus integrated spectral classification 

We used the standard optical diagnostic diagrams (BPT dia- 
grams, Baldwin, Phillips, & Terlevich 1981) to classify galax- 
ies into the H n-like, LINER, and Seyfert spectral types. For the 
spectral classification of the nuclear and integrated activity we 
used two different sets of boundaries in the optical line ratio di- 
agrams. The first are the classical semi-empirical boundaries of 
Veilleux & Osterbrock (1987, V&087 hereafter), which have the 
added advantage of making the comparison with previous results 
on these LIRGs easier. The V&087 boundaries are shown in 
Fig. 5. The second set includes the latest empirical and theoret- 
ical boundaries derived by Kauffmann et al. (2003) and Kewley 
et al. (2001a, b; 2006), and are shown in Fig. 6. Kewley et al. 
(2001a, b) modeled AGN and starburst line ratios to provide the- 
oretical boundaries in these diagnostic diagrams. They defined 
the so-called maximum starburst lines above which the line ra- 
tios cannot be explained by pure star formation. Kauffmann et al. 
(2003) and Kewley et al. (2006) derived empirical boundaries in 
these diagrams to separate H n, LINERs and Seyfert galaxies us- 
ing large samples of galaxies drawn from the Sloan Sky Digital 
Survey (SSDS). We will refer to this second set of boundaries as 
theoretical/SSDS. 

We did not attempt to correct the line ratios for extinction. 
The effects of extinction on the spectral classifications should be 
moderately small because the lines involved in the line ratios are 
in most cases close in wavelength. 

The first result worth noticing is that when using the 
[Om]^5007 fH/3 vs. [Nn]/f6584/Ha diagram and the theoreti- 
cal/SSDS boundaries (Fig. 6 and Table 3), a large fraction of 
the LIRG nuclei are classified as composite, and thus they are 
likely to contain a metal-rich stellar population and an AGN 
(see Kewley et al. 2006). It is also possible that some of these 
composite objects have an added contribution from shock ex- 
cited emission (from supernovae) associated with an aging star- 
burst (see Alonso-Herrero et al. 2000b). If we use all three 



BPT diagrams the fraction of composite objects appears to 
be smaller. However, Kauffmann et al. (2003) did not pro- 
vide an empirical separation between AGN and star-forming 
galaxies in the diagrams involving the [Oi]/I6300/Ha and the 
[S n]AA61 17,673 1/Haratios, so we cannot assess whether galax- 
ies are composite or not using those two diagrams. 

Using all three BPT diagrams and the V&087 classical 
boundaries (thin lines in Fig. 5) we find that only two nuclei 
in our sample of LIRGs show pure H n-like excitation (that is, 
classified as H n using all three BPT diagrams), five have an in- 
termediate LINER/H n classification, two fall in the LINER re- 
gions, and one nucleus is classified as a Seyfert galaxy (upper 
panels of Fig. 5, see also Table 3). This large fraction of com- 
posite objects (i.e., star formation and AGN activity) is a well- 
known property of samples of infrared-bright galaxies (Veilleux 
et al. 1995; Kewley et al. 2001b; Chen et al. 2009), regardless of 
what boundaries are used for classifying galaxies. 

We cannot provide a conclusive classification for IC 860 be- 
cause the emission lines in the blue part of the spectrum are 
not detected (note the strong H/3 absorption, Fig. 4). The high 
[N n]A65U/Ha and [S n] ^67 17, 67 3 1/Ha line ratios could indi- 
cate a Seyfert or LINER classification, but the correction for un- 
derlying stellar Ha absorption would decrease significantly the 
observed line ratios in this galaxy. 

A large fraction of the LIRGs in this sample show a strong 
contribution from an evolved stellar population as indicated by 
the presence of absorption features (Balmer line series, H+K 
Call lines) in the blue part of the spectra (Fig. 4). Although a 
detailed modeling of the stellar populations is beyond the scope 
of this paper, we can attempt to correct the [O m]/i5007/Hj8 ratio 
for the presence of an evolved stellar population for those galax- 
ies located near the AGN/HII boundaries. These are NGC 23, 
UGC 1845, and NGC 7591 (Fig. 6). Our preliminary model- 
ing indicates that a combination of intermediate (1-5 Gyr) and 
young (< lOMyr) stellar populations with a measured equiva- 
lent width of H/3 in absorption of ~ 4A would produce accept- 
able fits to their spectra. This is within the average H/3 stellar ab- 
sorption corrections obtained by Moustakas & Kennicutt (2006) 
for a sample of nearby star-forming galaxies. The approximate 
effects of correcting the observed [Oin]/l5007/H/3 for underly- 
ing stellar absorption are shown as arrows in the upper panel 
of Fig. 6. The corrections for the other line ratios involving Ha 
would be smaller. As can be seen from this figure, this correc- 
tion does not change fundamentally the result that these nuclei 
appear to be composite in nature. 
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Table 3. Spectral classifications of the nuclear and integrated spectra. 



Galaxy 


Type 


Classical V&087 Boundaries 


Theoretical/SSDS Boundaries 


Adopted 






[01] 


[Nil] 


[sn] 


[01] 


[Nil] 


[SII] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 23 


Nuclear 


L 


L 


L 


HII/L 


AGN/Composite 


HII/L 


Composite 




Integrated 


HII 


L/HII 


HII 


HII 


Composite 


HII 


HII 


MCG +17-07-001 




HII 


HII 


HII 


HII 


Composite 


HII 


HII 




Integrated 


HII 


HII 


HII 


HII 


Composite 


HII 


HII 


UGC 184*) 




L 


L 


L/HII 


Sy/L 


AGN 


Sy/HII 


Cnmnn^itp 




Integrated 


L/HII 


L 


HII 


Sy/HII 


AnN/Pnmnn';] tp 


HII 


Composite 


NGC 2388 


Nuclear 


HII 


L/HII 


HII 


HII 


Composite 


HII 


HII 




Integrated 


HII 


L/HII 


HII 


HII 


Cnmnrwitp 


HII 


HII 


MCG +02-20-003 


Nuclear 


L/HII 


L/HII 


L/HII 


HII 


Composite 


HII 


HII 




Integrated 


HII 


HII 


L/HII 


HII 




HII 


HII 


IC860 


Nuclear 




L/Sy? 


L/Sy? 




AGN? 


L? 




NGC 5939 


Nuclpar 


HII 


L/HII 


HII 


HII 


Composite 


HII 


HII 




Integrated 


HII 


L/HII 


HII 


HII 


HII 


HII 


HII 


NGC 6701 


Nuclear 


L 


L 


L/HII 


HII/L 


Composite 


HII 


Composite 




Integrated 


L/HII 


L 


L/HII 


HII 


Composite 


HII 


HII 


NGC 7469 


Nuclear 


Sy 


Sy 


Sy 


Sy 


AGN 


Sy 


Sy 




Integrated 


Sy/HII 


Sy 


Sy/HII 


Sy/HII 


Composite 


Sy/HII 


Composite 


NGC 7591 


Nuclear 


L 


L 


L 


L 


AGN/Composite 


HII/L 


Composite 




Integrated 




L/HII 


L/HII 




AGN/Composite 


L/HII 


Composite 


NGC 7771 


Nuclear 


HII 


HII 


HII 


HII 


HII 


HII 


HII 




Integrated 




L/HII 


L/HII 




Composite 


HII 


HII 



Notes.— Column (2): PMAS type of ID spectra. Columns (3), (4), and (5): Spectral classification based on the [O m],l5007/H/J vs. [O i]i6300/Ha, 
[Nn]/i6584/Ha', and the [S n]/U67 17,673 1/Ho- diagrams, respectively and the classical boundaries of V&087 in the BPT diagrams. The clas- 
sifications are, HII=HII galaxy, L=LINER, and Sy=Seyfert. Galaxies classified as Composite are likely to be a combination of AGN activity 
and star formation. The classifications take into account the ±l<x uncertainties reported by Kewley et al. (2001a) for the maximum starburst 
lines. Columns (6), (7), and (8): As Columns (3), (4) and (5) but using the latest theoretical/SSDS boundaries. Column (9): Adopted spectral 
classification. 



We have seven galaxies in common with the work of Veilleux 
et al. (1995). Our nuclear classifications are in good agreement 
with theirs. The only two exceptions are NGC 23 and NGC 6701, 
which we would classify as LINER/HII with the V&087 bound- 
aries, or as composite using the theoretical/SSDS boundaries. As 
explained in Sect. 3.2 we extracted our nuclear spectra with the 
smallest possible physical sizes allowed by the PMAS spaxels 
(~ 300pc, Table 2), whereas Veilleux et al. (1995) used linear 
sizes of 2 kpc for extracting their nuclear spectra. As can be seen 
from Fig. la and lh, the Veilleux et al. (1995) apertures included 
a large number of H n regions in the ring of star formation of 
NGC 23 and in the inner spiral structure of NGC 6701. This 
readily explains the H n-like classification given by Veilleux et 
al. (1995). 

The BPT diagrams for the integrated emission over the 
PMAS FoV are presented in the lower panels of Figs. 5 and 
6. The integrated line ratios of the four galaxies whose nu- 
clei are classified as LINERs using the V&087 boundaries now 
fall in the H n region or in the intermediate LINER/H n region. 
A similar situation is seen for NGC 7469 for which the inte- 
grated line ratios move toward the composite area in all dia- 
grams. This is well understood in terms of the increased con- 



tribution of extra-nuclear high surface-brightness Hn regions 
(see Ha morphologies in Fig. 1) to their integrated emission. 
For the other galaxies there is no general trend, as the inte- 
grated line ratios depend on the relative contribution of H n re- 
gions and diffuse emission to the total line emission over the 
PMAS FoV (Alonso-Herrero et al. 2009, in preparation). For 
instance, the integrated line ratios of NGC 7771 show larger 
[Nn]/16584/Ha and [S n]/l/l6717,6731/Ha line ratios than the 
nuclear ones, whereas the nuclear and integrated line ratios of 
other galaxies (e.g., MCG+02-20-003) remain approximately 
constant or slightly more similar to those of H n regions. 

Finally, we show an example of the power of optical IFS 
in identifying different excitation conditions. For NGC 7469 we 
also extracted the spectrum of an Hn region in the circumnu- 
clear ring of star formation, located at about 2" west from the 
nucleus (Fig. 4). It is clear that the line ratios are typical of H n- 
like excitation and are not contaminated by the nearby Seyfert 1 
nucleus (i.e., no broad components are present in the hydrogen 
recombination lines). 

Summarizing, the comparison of the nuclear and integrated 
activity classifications, together with the spatial distribution of 
the bright emission lines (in particular Ha) allowed us to isolate 
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Table 4. Gas kinematic results. 



oaiaxy 




Al! 


^*nuc,Ho' 


(1) 


(A) 


(J) 


(A\ 

w 


NGC 23 


4495 ± 9 


-215, +236 


139 + 17 


MCG +12-02-001 


4732 ± 8 


-110, +156 


105 + 14 


UGC 1845 


4796 ± 9 


-237, +181 


142 ± 17 


JN(j(_ Z3oo 


4U93 ± 9 


— 13 /, + 163 


1 ir . if 

1 13 ± 13 


MCG +02-20-003 


4950 ± 9 


-141, +110 


96 + 13 


IC 860* 


4146 ± 13 






NGC 5936 


3984 ± 13 


-66, +179 


66 ± 15 


NGC 6701 


3946 ± 13 


-41, +102 


92 ± 16 


NGC 7469 


4946 ± 9" 


-110,+207" 


125 ± 16" 


NGC 7591 


4869 ± 9 


-179, +345 


144 ± 17 


NGC 7771 


4334 ± 9 


-181, +199 


106 ± 14 



Notes. — Column (2): cz n uc,Ha is the Ha velocity of the nucleus, defined 
as the peak of the 6200 A continuum, except for NGC 7771 where it is 
refered to the position that makes the velocity field gradient symmetric 
(see text and Fig. 3c). Column (3): Avho- is the Ho- peak-to-peak velocity 
range over the PMAS FoV. Column (4): <Tn UC ,H a is the nuclear Ha ve- 
locity dispersion (corrected for instrumental effects). All measurements 
are in km s -1 , and are not corrected for inclination and/or angular reso- 
lution effects. 

*For IC 860 cz nuc is from [Nn]/16584 and [Sn],U6717,6731. We note 
that NED quotes a value of 3347 km s _1 from 21 cm H i measurements, 
whereas our value of the velocity is confirmed by measurements of the 
mid-infrared emission lines with Spitzer/IRS (Pereira-Santaella et al. 
2009, in preparation). 

"For NGC 7469 the measurements are from [Nn]/16584. 



and quantify the different ionization sources (nuclei, and circum- 
nuclear H n regions and diffuse emission) contributing to the ob- 
served emission in galaxies (see also Alonso-Herrero et al. 2009, 
in preparation). 

6. Ionized Gas Kinematics 

The velocity fields of the emission lines over the central few kpc 
(typically 5 kpc) are mostly consistent with rotation (see Fig. 3) 
for all the galaxies in the sample. We note that even on the phys- 
ical scales probed by PMAS (~ 300 pc), some velocity fields 
appear to show some peculiarities. However, it is clear that these 
velocity fields are more similar to those of disk galaxies (Falcon- 
Barroso et al. 2006: Daigle et al. 2006) than to those of local 
ULIRGs (Colina et al. 2005: Monreal-Ibero et al. 2006). The 
large scale continuum major photometric axis and the Ha major 
kinematic axis are broadly in agreement with each other (seen 
in projection) in the majority of our LIRGs. The ionized gas 
peak-to-peak velocities of LIRGs (over the central ~ 5 kpc) are 
typically between 200 and 400 km s _1 (see Table 4), whereas in 
ULIRGs with tidally induced flows the gas peak-to-peak veloci- 
ties can be as high as 600kms _1 (Colina et al. 2005). 

Three galaxies in our sample show circumnuclear rings of 
star formation. The central Ha velocity field of NGC 23 is con- 
sistent with rotation, with the major kinematic axis aligned with 
the major axis of the galaxy rather than with the orientation of 
the ring of star formation as seen in Ho- (Fig. la). A similar situ- 
ation is observed for NGC 7469, where the [Nn]/16584 velocity 
field has a major axis similar to that of the continuum, although it 
is asymmetric towards the north-west direction. This asymmetry 
was already reported in the rotation curve along the major axis 
of the galaxy by Marquez & Moles (1994). In both galaxies, it is 
likely that most of the mass in the central region is in relatively 
evolved stars rather than in young ionizing stars (see also Diaz- 



Santos et al. 2007 for NGC 7469). The Ha velocity field of the 
inner ~ 10" of NGC 7771 has a major kinematic axis aligned 
with the approximate orientation of the circumnuclear ring of 
star formation (east-west direction), whereas the outer velocity 
field appears to have a major kinematic axis in better agreement 
with that of the large scale continuum emission. This kind of ve- 
locity fields with symmetric distortions can be associated with 
the presence of a warped disk. It is also worth mentioning that 
there is dynamical evidence that NGC 7771 is weakly interact- 
ing with NGC 7770 (Keel 1993) and is located in a group of 
galaxies. 

Another example of a central Ha velocity field that deviates 
from perfect rotation is that of NGC 6701. This galaxy shows a 
complex overall morphology, with the presence of a large scale 
bar, an inner isophote twist produced by a spiral like inner ring 
(see Marquez et al. 1996, and the continuum 1.6/im image in 
Fig. lh), a perturbed rotation curve, and a small companion 
likely to be responsible for some of these properties (Marquez 
et al. 1996). As seen in other galaxies in our sample, the ma- 
jor kinematic axis of the Ha velocity field appears to be better 
aligned with the photometric axis of the galaxy (PA= 1 14°, Vogt 
et al. 2004) rather than with the orientation of the nuclear high 
surface-brightness emission Paa emission, which is seen in an 
almost north-south orientation. 

Arp 299 the most luminous galaxy in our sample and a 
strongly interacting system shows very complicated velocity 
fields (both of neutral and ionized gas, see Garcfa-Marin et al. 
2006), not only at the interface between the two galaxies, but 
also in the nuclear regions of the two members of the system. 
Moreover, the velocity fields of the ionized gas in Arp 299 do 
not appear to be dominated by ordered virialized motions, prob- 
ably as a result of the interaction between the two galaxies, as is 
the case of most ULIRGs. 

As can be seen from Fig. 3 for most LIRGs in our sample the 
peak of the Ha velocity dispersion coincides with the peak of the 
optical and near-infrared continuum emission (i.e., the nucleus), 
and thus the velocity dispersion is likely to be tracing mass. The 
nuclear Ha velocity dispersions are between 66 and 144km s _1 
(Table 4). The largest velocity dispersions of the ionized gas are 
associated with three of the nuclei classified as LINERs and the 
Seyfert galaxy NGC 7469. 

In the case of galaxies with circumnuclear rings of star for- 
mation, we find that the gas velocity dispersions in the rings are 
less than in the nuclear regions. This is consistent with the find- 
ings of Falcon-Barroso et al. (2006) for the same type of galax- 
ies and was interpreted as an indication for the presence of large 
amounts of cold gas from which stars have recently formed. 

There are only a few measurements of the stellar velocity 
dispersion of LIRGs. From near-infrared CO absorption features 
the stellar velocity dispersions are between 60 and 160kms _1 
(Shier, Rieke, & Rieke 1996; Hinz & Rieke 2006), similar to the 
range of Ha velocity dispersions measured for our LIRGs. The 
only LIRG in our sample with a stellar velocity dispersion mea- 
surement is NGC 7469, for which Onken et al. (2004) obtained 
cr* = 142 + 16kms _1 from the calcium triplet, in good agree- 
ment with our measurement from the ionized gas. The velocity 
dispersions of LIRGs are in general similar to or slightly less 
than the typical nuclear gas and stellar velocity dispersions of 
ULIRGs (Colina et al. 2005; Dasyra et al. 2006). This suggests 
that at least some LIRGs have comparable dynamical masses to 
those of ULIRGs, although a detailed modeling is needed to as- 
sess this issue. 



12 



Alonso-Herrero et al.: PMAS integral field spectroscopy of local LIRGs 



7. Discussion and Summary 

This is the first paper in a series presenting PMAS optical IFS 
observations of the northern hemisphere portion of the volume- 
limited (v = 2750 - 5200 km s _1 ) sample of local LIRGs defined 
by Alonso-Herrero et al. (2006). This sample is in turn part of 
the larger IFS survey of nearby (z < 0.26) LIRGs and ULIRGs 
assembled by Arribas et al. (2008). In this paper we presented 
the observations and data reduction of the PMAS observations. 
The PMAS observations cover the central 16" x 16" (typically 
the central 5 kpc) with spaxels of 1" in size, and a spectral range 
~ 3800 - 7200 A. The PMAS IFS data were complemented with 
our own existing near-infrared HSr/NICMOS observations of 
the 1.6 fim continuum and the Paa emission line. The main goal 
of this paper is to present an atlas of the observations and the 
general IFS results of the sample of LIRGs and compare them 
with local ULIRGs. 

On the physical scales probed by the PMAS IFS (~ 300 pc) 
the optical and near-infrared stellar morphologies are similar for 
most galaxies, indicating that extinction is not playing a major 
role for this sample of LIRGs, except in the innermost regions. 
Similarly, there are no major morphological differences between 
Paa and Ha. The HST/NICMOS Paa and PMAS Ha observa- 
tions are complementary, with the former revealing in great de- 
tail (physical scales of a few tens of parsecs) the morphologies 
of the high surface-brightness H n regions, and the latter being 
sensitive not only to H n regions but also to diffuse emission with 
lower surface brightness. 

In the majority of the LIRGs in our sample the peaks of the 
continuum and gas (e.g., Ha, [Nn]/16584) emissions coincide. 
This contrasts with local interacting ULIRGs, where the extreme 
effects of the interaction processes on the ionizing mechanisms 
and the distribution of dust can cause displacements between 
the peaks of continuum and gas emission of typically 2-4 kpc 
(Colina et al. 2005: Garcfa-Marin et al. 2009). The only excep- 
tions in the LIRG sample are galaxies with circumnuclear rings 
of star formation (NGC 23 and NGC 7771) and the strongly in- 
teracting galaxy Arp 299. In the case of the galaxies with cir- 
cumnuclear rings of star formation, the most luminous Ha emit- 
ting regions are found in the rings rather than in the nuclei of the 
galaxies, and the displacements are well understood in terms of 
differences in the stellar populations. In Arp 299 the displace- 
ments are due to the high extinctions suffered by the nuclear 
regions (see Alonso-Herrero et al. 2000a and Garcfa-Mann et al. 
2006). 

Using standard BPT diagrams we compared the excitation 
conditions of the nuclear and integrated (over the PMAS FoV) 
ID spectra of the LIRGs in our sample. Only two nuclei show 
pure H n-like excitation using the classical V&087 boundaries, 
and one has a Seyfert nucleus. The rest are classified as LINER 
or intermediate LINER/Hn, using the V&087 boundaries, or 
alternatively as composite objects using the theoretical/SSDS 
boundaries. We also found that a large fraction of the nuclei 
show evidence for a strong contribution from an evolved stel- 
lar population (absorption features). There is no general trend 
for the excitation conditions of the integrated emission when 
compared with the nuclear excitation, as the former depends 
on the relative contributions of Hn regions and the diffuse 
emission to the line emission over the PMAS FoV (Alonso- 
Herrero et al. 2009, in preparation). That is, galaxies dominated 
by high surface-brightness Hn regions show integrated Hn- 
like excitation, whereas galaxies with more diffuse low surface- 
brightness emission tend to show slightly larger [Nn]/16584/Ha 
and [S ii]/L16717,6731/Ha line ratios. 



The Ha velocity fields over the central few kpc covered 
by the PMAS observations are generally consistent, at least to 
first order, with rotational motions. The observed Ha velocity 
amplitudes (peak-to-peak) are between 200 and ~ 400kms~ 1 . 
Although the velocity fields of some LIRGs show some peculiar- 
ities, they are not as perturbed as those of most local, strongly in- 
teracting ULIRGs. In that respect, the velocity fields of the emis- 
sion lines of our sample resemble those of disk galaxies (Falcon- 
Barroso et al. 2006: Daigle et al. 2006). In most LIRGs in our 
sample the peak of the Ha velocity dispersion coincides with the 
peak of the optical and near-infrared continuum emission (i.e., 
the nucleus). Thus, the velocity dispersion is likely to be tracing 
mass and provides further support for rotation. The nuclear Ha 
velocity dispersions are in the range <x nucHff = 66 - 144kms~ 1 
and are similar to the stellar values measured from near-infrared 
CO absorption features measured for other LIRGs. The LIRG 
nuclei with the largest Ha velocity dispersions are those classi- 
fied as LINERs and the Seyfert 1 nucleus of NGC 7469. 

Throughout this paper we discussed the ionized gas and stel- 
lar distributions, excitation conditions and kinematics of a vol- 
ume limited sample of LIRGs. We also showed that the proper- 
ties of this sample of LIRGs, and in particular, their kinematics, 
are more similar to those of disk galaxies, rather than to those 
of local, strongly interacting ULIRGs. Our LIRGs are part of a 
flux and volume limited sample drawn from the IRAS RBGS of 
Sanders et al. (2003). As a consequence, our sample is predom- 
inantly composed of low-luminosity LIRGs with an average IR 
luminosity of log(LiR/L ) ~ 11.32 for the full sample (north- 
ern and southern hemispheres, see Alonso-Herrero et al. 2006). 
Sanders & Ishida (2004) demonstrated that log(L IR /L ) ~ 1 1.50 
marks the transition between samples being dominated by disk 
galaxies and merger dominated samples. In the Sanders & Ishida 
(2004) sample of LIRGs, most objects at log(L IR /L ) < 11.40 
are spiral galaxies with no signatures of a major interaction and 
pairs of galaxies, whereas at log(Li R /L ) > 11.70 most objects 
are strongly interacting equal mass galaxies with overlapping 
disks. The similarities of our sample with normal disk galax- 
ies are then well understood because our sample is mostly com- 
posed of disk galaxies and galaxy pairs not undergoing a strong 
interaction (Sect. 2.1). 

At intermediate redshifts (z ~ 1) LIRGs are the main con- 
tributors to the star formation rate density (Elbaz et al. 2002; 
Le Floc'h et al. 2005; Perez-Gonzalez et al. 2005; Caputi et al. 
2007). Moreover, these z ~ 1 LIRGs are mostly classified as spi- 
ral galaxies (Bell et al. 2005: Melbourne et al. 2005, 2008) with 
the star formation smoothly distributed in the disks of the galax- 
ies as seen for our sample of LIRGs. There is also a strong evolu- 
tion in galaxy kinematics at z < 1, but still about half of the line 
emitting population at these redshifts, which are mostly LIRGs, 
have rotating disks (Puech et al. 2008; Yan et al. 2008). While 
local and distant LIRGs may have self regulated star formation 
as in disk galaxies (see Bell et al. 2005), it is clear that their star 
formation rates are at least a factor of ten or more higher than in 
spirals. Locally the central regions of LIRGs are found to contain 
a large population of giant HII regions not observed in normal 
spiral galaxies (see Alonso-Herrero et al. 2006). This, together 
with a higher star formation efficiency of the dense gas in local 
LIRGs (Gracia- Carpio et al. 2008) may explain the high star for- 
mation rates of LIRGs when compared to normal disk galaxies. 
Studying the spatially resolved properties of flux-limited com- 
plete samples of local LIRGs may help us understand if the sim- 
ilarities with z ^ 1 LIRGs stem from the same physical processes 
and comparable evolutionary states. 
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Summarizing, we demonstrated the ability of optical IFS to 
spatially resolve the different ionization sources contributing to 
the observed emission of LIRGs, as well as to study their kine- 
matic properties. Full detailed studies of the extinction, exci- 
tation conditions, stellar populations, and kinematics of LIRGs 
and ULIRGs will be presented in forthcoming papers. 
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Fig. 4. PMAS spectra (plotted in arbitrary units) of the nuclei of the galaxies in the sample. Additionally, for each LIRG the inset 
shows the blue part of the integrated spectrum so that the absorption features as well as the H/3 and the [Om]/15007 emission lines 
can be clearly seen. In the case of NGC 7469 we also show in the lower right panel the full spectrum of an H n region located in the 
circunmnuclear ring of star formation. 
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Fig. 4. Continued. 
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Fig. 5. BPT diagrams for the PMAS nuclear (upper panels) and the integrated (lower panels) emission of the LIRGs in our sample. 
The thin lines are the classical boundaries of Veilleux & Osterbrock (1987) for the H n region, LINER, and Seyfert excitation. 
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Fig. 1. (a) NGC 23. The middle and bottom panels are the PMAS observed (not corrected for extinction) maps of the brightest 
emission lines: H/3, [Om]/t5007, [Oi]^6300, Ha, [Nn]A6584, and the sum of the [S n]^6717, 6731 lines. The H/3 map has not 
been corrected for stellar absorption. The upper left panel is the map of the ffiiT/NICMOS Paa emission line. The maps of the 
PMAS 6200 A and the //ST/NICMOS 1 .6 fim continuum emission are the upper right and the upper middle panels, respectively, 
both representing the stellar emission. The cross on the PMAS maps shows the location of the PMAS 6200 A continuum peak. Due 
to the lack of absolute astrometry of the PMAS observations, the PMAS continuum peak is not shown on the NICMOS images. The 
horizontal bar in the Pao- panel represents the 2kpc linear scale used by Kim et al. (1995) and Veilleux et al. (1995) for extracting 
their nuclear spectra. The contour shown on the PMAS continuum map corresponds to the external isophote used for extracting the 
integrated ID spectra (see Sect. 3.2). The orientation of the images is north up, east to the left. All the images are shown in a square 
root scale. 



Alonso-Herrero et al.: PMAS integral field spectroscopy of local LIRGs 



19 



MCG+12-02-001 



1.6/xm I O 



o 

CD 

m O 

u 

cd 




Ik Lu 



5 0-5 5 

arcsec arcsec 
Fig. 1. (b) As Fig. la but for MCG +12-02-001. 
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Fig. 1. (f) As Fig. la but for IC 860. 



5 -5 5 -5 

arcsec arcsec 
The HjS, [Om]/15007 and [Oi]/I6300 emission lines are not detected in this galaxy. 
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Fig. 1. (i) As Fig. la but for NGC 7469; there is no HSTfNlCMOS NIC2 Pao- image available for this galaxy. The PMAS Wfi and 
Ha maps shown in this figure were constructed by fitting one component to the lines (see text and also Fig. 2). 
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Fig. 1. (j) As Fig. la but for NGC 7591. 
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Fig. 1. (k) As Fig. la but for NGC 7771. The PMAS mosaics were constructed with the east and west pointings done for this galaxy 
and they cover approximately the central 28" x 16" region. 



